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Abstract 



The application of the Gardner method for generation of conservation laws to 
all the ABS equations is considered. It is shown that all the necessary informa- 
tion for the application of the Gardner method, namely Backlund transforma- 
X/^ • tions and initial conservation laws, follow from the multidimensional consistency 

^ • of ABS equations. We also apply the Gardner method to an asymmetric equa- 

tion which is not included in the ABS classification. An analog of the Gardner 
method for generation of symmetries is developed and applied to discrete KdV. 
It can also be applied to all the other ABS equations. 

> 

(N '. 1 Introduction 

The first integrable example of a partial difference equation (PAE) goes back to 
Q ■ who wrote down an equation on a quad-graph related by a simple change of variables to 

O . the discrete KdV equation. However, the subject has only been intensively developed 

in the last few years. A particular milestone was the so-called ABS classification of 
integrable scalar quad-graph equations based on the principle of consistency on the 
^ ' cube [1] . All the equations in this classification have the discrete analogue of a matrix 

' Lax pair and also a natural auto-Backlund transformation [2]. Other integrability 

properties (discrete analogues of scalar Lax pairs, bihamiltonian structures and infinite 
numbers of conservation laws) have been studied less. 

The topic of conservation laws for the PAE has become very important recently. 
This started in the work of Orphanidis [12] in which he presents conservation laws for 
the fully discrete sine-Gordon equation. The next step in the research of conservation 
laws was the introduction of a systematic computation method. The first systematic 
method developed was the direct method [HI UHl [16]. This method allows use of 
computer algebra, which makes it more constructive. But it has some disadvantages: 
it requires massive computations and cannot produce infinite numbers of conservation 
laws. The next method for computing conservation laws was proposed in [17]. It 
produces infinite numbers of conservation laws by acting with a mastersymmmetry on 
basic conservation laws. 
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A new approach to the computation of conservation laws for PAE has appeared 
very recently [H]. It is called the Gardner method and is an analog of the Gardner 
method for partial differential equations. The method uses a Backlund transformation 
(BT) to generate conservation laws. It seems that Gardner method is the most efficient 
among the methods described above. 

Symmetries for PAE were also researched recently. They first appeared as sim- 
ilarity constraints for integrable lattices [lO]. Tongas et al. pointed out that the 
similarity constraints for quad-graph equations obtained previously are equivalent to 
characteristics of symmetries [20j. Symmetries of several quad-graph equations have 
been found in [TJ [SI [T31 HH]- Hydon developed a direct method for finding symmetries 
based on solving functional equations by creating an associated system of differential 
equations that can be solved O |6]. In [18] this method was applied to PAE. There 
the authors compute five-point symmetries for all the ABS equations (equations from 
PP). Symmetries for PAE were also discussed in other articles [HIEI]- 

Our intention is to show that the Gardner method for generation of conservation 
laws for PAE can be applied systematically. We apply it to all the ABS equations 
and to an asymmetric quad-graph equation. We also show that the Gardner method 
can be applied for symmetries. Namely, we use it to generate an infinite number of 
symmetries from known one. 



2 The Gardner method 

Before presenting the Gardner method for generating conservation laws for PAE, we 
review the method for the continuum Korteweg-de Vries equation (KdV). The Gardner 
method for KdV starts with the BT. The BT states that if u solves KdV then so does 
u + Vx where v is a. solution of the system 

Vx = 6 — 2u — 

Vt = -^Uxx + {0 + U) [9 - 2u - V^) + UxV 

Here ^ is a parameter. It is straightforward to check that v is actually the G component 
of a conservation law. Specifically, we have 

dtv + dx (^Ux - (m + 0)v^ = . 

But V depends on the parameter 6. Expanding in a suitable series in 9 yields an infinite 
number of conservation laws. For KdV the appropriate series is 



il/2 ^ .'^x t^xx 2^ Uxxx 8uUx Uxxxx 8u -\- lOli^. + 12uUxx ^ 

29 4^3/2 ^2 16^5/2 + ^ ) ■ 



In [11] the implementation of the Gardner method to the dKdV equation was given. 
Here we apply the Gardner method to the following asymmetric quad-graph equa- 
tion [TT] 

tt^i0,0Wl,0 - /3(m0,0W0,1 + Wi,oMl,l) = 0. (1) 

Here k^l eT? are independent variables and mo,o = u{k, I) is a dependent variable that 
is defined on the domain 1?. We denote the values of this variable on other points by 
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Uij = u{k + + j) = SlSjuofi, where Sk, Si are the unit forward shift operators in k 
and / respectively. Equation ([1]) is not included in the ABS list, since it is asymmetric. 
A BT for is 

a{uo,oUi,o + UofiUifi) - 0{uQfiUofi + ni^o^ti.o) = 0, (2) 
OuofiUQfi - /3(mo,oMo,i + uoflUo,i) = 0. (3) 

Here ^ is a parameter. This BT follows from the multidimensional consistency of 
equations (fT|2|3|) embedded in three dimensions. Let us introduce the third variable 
m so that 

Uk,i,m = %,o,o = %,o, ^io,o,i = %,o, Mi,o,i = Mi,o, etc. (4) 
Then, equations ffTP|5]) can be written as 



Pi 



awo,o,oMi,o,o - /3(mo,o,oMo,i,o + mi,o,o«i,i,o) = 0, 

a(%,o,oMi,o,o + Wo,o,i^i,o,i) - 6'(mo,o,oWo,o,i + "1,0,0^^1,0,1) = 0' (5) 
6*^^0,0,0^^0,0,1 - P{uofl,oUo,ifl + ^0,0,1^^0,1,1) = 0. 



Equations P2 and P3 form a BT for Pi. In the sequel we use the notation (|H), writing 
Mo,o,i and wo,o,-i for the BT and inverse BT of ^0,0,0, where possible, however, we revert 
to 2 dimensional notation. 

In general we cannot solve the equations of the BT to write ^0,0,1 iii terms of tto,o,o- 
However there is a special case ^ = a for which mo,o,i = 'iii,o,o or m_i,o,o- Consider 
6 = a + e where e is small, and look for a solution of BT in the form 



00 



uofi,i = uifi + ^ v^Qle\ (6) 

i=l 

We just look at the first equation of the BT. This reads 

00 / CXD \ / 00 

i=i V j=i / \ i=i 

The leading order approximation gives 

^(1) ^ "l,o("0,0 + U2,o) 

a(w2,o - Uo,o) 

Higher order terms give 

(i-l) , (i-l) sr^i-l (i) (j-i) 

Voo = 7 H , t = 2,3.... (7) 

"("2,0 - Wo,o) 

As in the case of the dKdV equation all these formulas are on the horizontal line, i.e. all 
the Vqq only depend on values of Uij with j = 0. An infinite sequence of conservation 
laws can be obtained starting from the e expansion of a single conservation law. It is 
straightforward to check that if we define 

F = lnf^V G = -lnf ""°'°'^~^"^'°'° V (8) 
V«o,o,oy V «o,o,o J 
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then 

iSk-I)F+{Si-I)G = 0, 

on solutions of ([5]). By plugging and 6 = a+e into ([8]) and expanding F = J2ilo 
and G = — ln(2e) + J2ilo Gi^^ we obtain 



Uo,0/ ^2,0-^0,0 

Go = - In ^1'° ^ - "i,o(3w2,o + uofi) + M3,o(m2,o - Wo,o) 



4n2,o%,o(%,o + ^^i,o) + (^1,0 - "o,o)(%,o - ^i,o) 

^ 2 

2(^2,0 - liO,o)^(«l,0 - ^i3,o) 

^ _ (^3,0^2,0 + 3m2,0Mi,0 + Mi.qMq.o - ^3,0^0,o)^ 

8(m2,0 - «0,o)^(^i3,0 - ^^l,o)^ 

Thus we see how the expansion of the BT around the point 6 = a yields an infinite 
sequence of conservation laws on the horizontal line. Expansion around 9 = j3 does 
not seem to yield an infinite sequence of conservation laws on the vertical line, since 
(II]) is asymmetric. In the case of the dKdV equation expansion around 6 = (3 also 
gives an infinite number of conservation laws. 



3 The Gardner method for conservation laws of in- 
tegrable equations on the quad-graph 

In this section we present all the necessary information for the application of the 
Gardner method to the quad-graphs that are listed in pQ. The general form of ABS 
equations is 

Pi : P(Mo,o,o,wi,o,o,'"o,i,o,^^i,i,o,a,/?) = 0. (9) 

For the application of the Gardner method we need a special BT and the initial 
conservation law. For all the ABS equations the necessary BTs are known and these 
are so called natural auto-BTs The general form of the natural auto-BTs for ABS 
equations can be obtained from the form of the equation For dH]) the form of the 
natural auto-BT is 

P2 ■ P{uoflfl,uiflfi,uofl^i,ui^o,i,a,6) = 0, (10) 
P3 ■ -P(%,o,0) %,o,i' ^0,1,0, «o,i,i' ^) = 0- (11) 

This is the BT which we need for the application of the Gardner method. Note that 
equations flU|10lllip are one equation embedded in three dimensions. 

The initial conservation law(ICL) also follows from multidimensional consistency. 
Expressions for conservation laws for Pi, P2, P3 are respectively 

{Sk - + {Si - I)Gi = 0, 

(Sk - I)F2 + (S^ - I)H2 = 0, (12) 
{Si - I)G, + (S^ - I)Hs = 0. 
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Five-point conservation laws for ABS equations were found in [T7], where the authors 
show that each ABS equation has three five-point conservation laws. Two five-point 
conservation laws for Pi can be always presented as 

F = /(% -1,0, Mo,o,o> '"0,1,0, G = g{uo-ifl, wo,o,o, ''^1,0,0, «, (3), (13) 

F' = 5f(M_i,o,o, tto,o,o, ^^0,1,0, 13, a), G' = /(m-i,o,o, ^^o,o,o, "1,0,0, «)• (14) 

As we said before, equations Pi, P2, P3 are one equation which is embedded in three 
dimensions. Therefore conservation laws for Pi, P2 can be obtained from conservation 
laws for Pi by rewriting them on the appropriate plane. The conservation law for P2 
which corresponds to F, G is 

F2 = /(mo,o,-i,Wo,o,o, "0,0,1,^), H2 = fi'(Mo,o,-i,%,o,o,^^i,o,o,a,6'). (15) 
The conservation law for P3 which corresponds to F, G is 

G3 = /(mo,o,-i, Wo,o,o, Wo,o,i,^), H3 = 5'(mo,o,-i,Mo,o,o,%,i,o,/5, 6*). (16) 

P2 is equal to G^, so with the help of linear combination of expressions from f ll2p we 
obtain that 

{Si - I){Sk - /)P2 + {Si - I){Sm - I)H2 - {Sk - I){Si - I)Gs 

- {Sk - I){Sm - I)Hs = {Srn - I){{Si - I)H2 - {Sk - I)Hs) = 

is true on solutions of P2 and P3. This expression can be integrated with respect to m 
to give 

{Si - I)H2 - {Sk - I)Hs = C{k, I) (17) 

where C{k,l) is a function which has to be determined. Equations P2 and P3 are 
symmetric with respect to the reflection of the coordinate m, that is 

P{UOfl,Q, ^1,0,0, ''^0,0,1, '«1,0,1, «, 9) = ±-P(%,0,l, "^^1,0,1, "^^0,0,0, '«1,0,0, «, 

= ±SmP{uofi,o, ^ii,o,o, Mo,o,-i, ^ii,o,-i, a, 0). 
Therefore ( !T7|) is also symmetric with respect to the reflection of coordinate m, so 

{Si - I)9{uo,o,i,Uo,o,o, Mi,o,o, a, ^) - {Sk - ^)5'(%,o,i, ^0,0,0, "0,1,0, d) = G{k, I). (18) 

By substitution of P2 and P3 into (fT8|) we obtain that C{k,l) = for all the ABS 
equations. Therefore 

{Si - I)H2 - {Sk - I)H, = 0, 

for all the ABS equations. Thus, the components of ICLs for all the ABS equations 
are 

FicL = — fi'(wo,o,i, ""0,0,0, ""0,1,0, tt, 0)1 GicL = fi'(«o,o,i, ^0,0,0, ""1,0,0, /3, 0). 

We present a table with ICLs for all the ABS equations bellow. The equations from 
the ABS classiflcation are as follows; for convenience, we have used the form of Q4 
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that was discovered by Hietarinta [3]. 



Ql : Q!(uo,o - Mo,i)(^i,o - "la) - Piuofl ~ ui,o)(woa - + S'^al3{a - /3) = 0, 

Q2 : a{uo,Q ~ uoa)('«i,o - "1,1) " ^(""o.o - ui,o)("o,i - ""1,1) 

+al3{a - ,3)(wo,o + Ui,o + Uo,i + Ui,i) - a/3(a - /3){a'^ - a/3 + /3^) = 0, 
Q3 : _ Q,2)(uo QUia + wi^oMo.i) + /^("^ - l)(uo,owi,o + uo,i"ia) 

-a(/32 - l)(uo,o"o,i + Mi.o^/i,i) - S^{a^ - f3^){c^^ - - l)/(4a/3) = 0, 

Q4 : sn(a)(uo, 0^1,0 + wo,i"ia) " sn(/3)(uo,oWoa + "1,0^1, 1) " sn(a - /3)(uo, 0^1,1 + wi^ouoa) 

+sn(a - ^)sn(Q!)sn(/3)(l + A'^Mo.oWi.oUoaUi.i) = Oj (19) 
HI : (uo,o - ui^i){ui^o - U04) + /3 - a = 0, 

H2 : (wo,o - ui,i)(ui,o - "0,1) + (/? - a)(wo,o + wi.o + mo,i + ""i,!) + - = 0, 

H3 : a(Mo,o'"i,o + wo,i"i,i) ^ ;^("o,o"o,i + ui^qUia) + S'^{a'^ ~ /3^) = 0, 

Al : a(uo,o + wo,i)(wi,o + "la) ^ P{uq,o + ui,o)(uoa + "1,1) - S'^al3{a - /3) = 0, 

A2 : (/32 - a2)(uo,oWi,o'«o,i"i4 + 1) + ^("^ - l)("o,oUoa + "i.c^i,!) 

-q;(/3^ - l)(wo,oUi,o + Mo,i"i,i) = 0- 

Here sn(a) = sn(a; K) is a Jacobi elliptic function with modulus K. Without loss 



of generality, the parameter 6 is restricted to the values and 1. Obtained ICLs are 
summarized in table 1, in which we list the components F and G for each of the ABS 
equations. 

Table 1: Initial conservation laws for the equations from the ABS classifica- 
tion 

Eq. Components 

F = In (-{to,o + uofl - 69) - \n{6{/3 -9) + - Uo,o) , 
G = - In (mo,o - uo,o - 69) + \n{6{a - 9) + wi,o - tto.o) , 

F = \n{{9p' + P{uo,o - uo,o - 9^) + 9{uo,o - ^io,i))') 
-ln(^^ - 2^2(^^0,1 + uofi) + (mo,o - Uq^i)^), 
G = - ln{{9a^ + a(Mo,o - Mo,o - d^) + OiuQ,o - Uifi))^) 
+ \ia{a'^ - 2a^{uifi + uo,o) + (mo,o - ^^1,0)^), 

F = ln((6'2(no,o - /3^^o,i))^ - ^(1 - a^)^to,o - a^Mo,o + a^Mo.i) 

- ln(4/3(/3wo,o - %,i)(%,o - + '^(l - 
G = - \n{{9'^{uofi - auifl) - 9uofi{l - a^) - a^Mo,o + aiti.o)^) 
+ ln(4a(awo,o - ^ii,o)(^^o,o - "^^1,0) + 6(1- a^)^), 

F = ln(sn(6')2(l + K'^uo,oUo,o) + 2cn{9)dn{9)uo,oUo,o - ^oo - ^00) 

- ln(sn(^ - /3)\uo,iUo,o - sn{9)snii3))isn{9)sn{f3)K^uo,oUo,i - 1) 
-(sn(6')Mo,i - sn(/3)zxo,o)(sn(/?)iio,i - sn(6')?io,o)), 
G = -ln(sn(^)2(l + K^uofiUofl) + 2cn(^)dn(^)no,o%,o - ""o.o - "^0,0) 
+ ln(sn(a — 6')^(mi,oWo,o — sn(6')sn(a))(sn(6')sn(a)A'^'Uo,o'iii,o — 1) 
-(sn(6')Mi,o - sn(a){to,o)(sn(a)iii,o - sn(6')Mo,o)), 



Continued on next page 
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Ql 
Q2 

Q3 
Q4 



Eq. Components 



HI 


F — 

G = 


illl"0,0 "0,1 J) 

ln({io,o - uifl), 


H2 


F = 
G = 


In {{l3 -9 + uo,i - Mo,o)^) - ln(/3 + uq^ + uq^i) , 
- In ((a - 6* + uifi - uoflY) + ln(a + uq^ + ui^) , 


H3 


F = 
G = 


In ((/3tto,o - 6'mo,i)^) - ln(5/3 + %,oWo,i) , 
- In ((a£to,o - Oui^oY) + ln(5Q; + %,o^ii,o) , 


Al 


F = 
G = 


ln((^io,o + ^o,o)' - S9^) - ln((i/o,i - ^o,o)' - ^{6 - /3)2), 
- ln((Mo,o + uo,o)^ - 66^) + ln((Mi,o - uo,o? - S{a - 9)^), 


A2 


F = 
G = 


ln({io,o%,o - G){GuQflUQfl - 1) - ln((^Mo,i - /3mo,o)(/3%,i - ^%,o)), 
- ln({to,oMo,o - 9){9uQflUQfl - 1) + ln((6'Mi,o - a{to,o)(awi,o - 6*^^0,0)), 



Another problem which can appear during the apphcation of the Gardner method 
is that Vq \ has to be found for i = 1,2, ... and for all the ABS equations. We prove 

that it is always possible to explicitly find \ for i = 1,2, ... for all the ABS equations 
in the application of the Gardner method in the k direction. For the / direction the 
proof is similar. Consider 6 = a + e where e is small, and look for a solution of BT in 
the form 

00 

^^0,0,1 = ^i,o + Xl^Se'- (20) 

1=1 

We look only at the first equation of BT. This reads 

(00 00 \ 

Mo,o, ^1,0, uifi + ^ v^'le\ U2fl + ^ v'fle', a,a + e] =0. (21) 
i=l i=l J 

Let Cn be the coefficient next to the e" in the expansion of (12T|) in the Taylor series 
around e = 0. At first sight C„ depends upon vl^j^, v^j^, i = 1,2, ...n and to find v^j^ one 
has to solve non-trivial difference equations. We show that that Cn depends just upon 
Vqq, Vil, i = 1,2, ...n, j = 1, 2, ...n — 1, so C„ can be solved explicitly with respect to 

Vq^J. First of all let us show 

Gq = P (mo,o, 1^1,0, Mi,o, M2,o, a, a) = 
Since P is an ABS equation it has the symmetry property [1] 

P{x, u, V, y, a, P) = -P{x, v, u, y, a, /3). 
According to this symmetry property we obtain 

P {uo,o, ui,o, Mi,o, U2,o, Oi, a) = -P (uo,o, ""1,0, ^^1,0, U2fl, a, a) = 0. 

The derivative of (1211) with respect e is 

,00 00 
^/ = I1 <^^'Ps + E ^vfle^-'P, + Pe, (22) 

1 = 1 4=1 
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where Pi is the derivative of P with respect to its ith. argument. For example 
d 

-P [ W0,0,Wl,0,«,«2,0 + 

i=l 



oo 



-P Uo,o, Mi,o, a, U2fl + ^ v^^le\ a,a + e 



a='"l,0+Ei=i ^0,0^ 



da 

For e = we obtain Ci 

'■ ^oflP^ (%,o, «i,o, Wi,o, ""2,0, a, a) + ^i^d A ("0,0, ""1,0, Wi,o, ^^2,0, «, «) 

+ -Pe (''^0,0, ^^1,0, Mi,o, M2,o, a, «) • (23) 



ae 



We checked that for all the ABS equations 

-P3 (^^0,0, Mi,o, Mi,o, W2,o, a, a) 7^ 0, P4 (mo,o, ^1,0, lii.o, ^^2,0, «, a) = 0. 

So, Ci depends upon Wq^q and does not depend upon w(^q. From ( l22l) it is seen that 
generally the form of Cn is similar to the form of Ci, namely 

Cn = nlv^'^QPs (mo,o, ^^1,0, uifl, U2fi, a, a) + nlv^^^^P^ {uofi, Mi,o, ui^, ^2,0, a,a) + R 

Here does not depend upon ijq^J, vj^J. As we said before P4 (mo,Oi 1*1,0, ^^i.Oj ""2,05 ct, a) = 

therefore C„ can be solved with respect to v^q. By solving Ci with respect to VqI 

iteratively for i = 1,2, ... we can find as many iJq^q as necessary. Once v^q, i = l...n 
are known, we plug (l20l) into the initial conservation law and expand it around e = 
up to order n. In this way n conservation laws can be found. 

4 The Gardner method for the symmetries of in- 
tegrable equations on the quad-graph 



In [18] five-point symmetries were found for ABS equations. The authors show there 
that each ABS equation Pi has four five-point symmetries. Two of these symmetries 
can always be presented as 

Xh = ?7(ii_i,o,o,wo,o,o, ""1,0,0, a)c?«o,o,o' (24) 
Xv = r]{uo-ifi,uofi^o,uo,ifl,l3)duo_o,o- (25) 

Therefore equations P2 and P3 both have a symmetry 

X = r]{uo,o-i,uofl,o,uofi^i,9)duo^o,o- (26) 

Since Pi, P2 and P3 are consistent in three dimensions we obtain that X is also a 
symmetry for Pi. Let us call X the initial symmetry for Pi. Let us apply the Gardner 
method to X or, in other words, expand X in a series. X involves the BT and inverse 
BT tto,o,i and -u^c-i, therefore we are looking for the solution of those for 6' = a ± e 
where e is small in the forms 

00 

%,o,i = «i,o,o + i^SJe*, (27) 
1=1 
00 



7/0,0,-1 = ^-1,0,0 + X^^oje*- (2^ 



i=l 
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By plugging this in the characteristic rj of the symmetry X we obtain 

oo oo 

r] = ?7(m_i,o,o + X] ^Se*, Wo,o,o> Wi,o,o + v^'^e', a + e). (29) 

i=l 1=1 

After expansion of r] in Taylor series around e = we obtain 

oo 

V = J2me'. (30) 

i=l 

Each r]i, i = 1,2,... is the characteristic of a symmetry for Pi. For example let us 
consider the discrete KdV equation 

(wo,o,o - wiA.o) (^1,0,0 - ^^o,l,o) + /3 -a = 0. (31) 
Two of the symmetries for this equation are 

Xh = f^«o,o,o' (32) 

"^1,0,0 ~ 

Xy = ^«o,o,o- (33) 

"^^0,1,0 " ^^0,-1,0 

Therefore 

X= (34) 

is also a symmetry. Expression ( l29l) for X is 

1 



We just look at the first equation of BT which is P2. This reads 



00 



E ^o!oe' j yuo,o,o - «2,o,o - ^S^y , (35) 

00 \ / 00 \ 

E "'Se^ ^-1,0,0 - ^1,0,0 + E . (36) 



.1=1 / \ i=l 

The leading order approximation gives 



<^ = , (37) 

%,o,o ~ ^2,0,0 

w[]l = ^ . (38) 

"^^-1,0,0 — "^1,0,0 



Higher order terms give 



« - ^ \ E-oUo'\ ^ = 2,3,... , (39) 



'1,0 — 

J=l 



i-1 

E^i^Ko'^ t = 2,3,.... (40) 



The first three coefficients in the expansion of rj around e = are 
1 



V3 



+ 



U-2,0 - "2,0 



+ 



U-3,0 - Uin 2 



(mO,0 - ^^-2,o)^(^i-3,0 - W-l,o)(«-l,0 - ^^l,o)^ (^iO,0 - W-2,o) (^^-2,0 " '"2,o)(W-l,0 " «l,o)^ 



These are characteristics of symmetries for ( 13T]) . rii and 772 were aheady presented in 

Eol. 



5 Concluding remarks 

In this article we showed how to apply the Gardner method for generation of conser- 
vation laws to all the ABS equation. It is shown that all the necessary information for 
the application of the Gardner method follows from the multidimensional consistency 
of ABS equations. Namely the natural auto-BT for ABS equation is the equation itself 
embedded in three dimensions. This construction is possible because of the multidi- 
mensional consistency of the equation. The initial conservation law for ABS equations 
follows from the linear combinations of conservation laws for the corresponding em- 
bedded equations. 

Another important result of this article is the introduction of the Gardner method 
for symmetry generation for ABS equations. In the Gardner method for the symmetry 
generation one has to substitute expressions 



Uo 



,0,1 = Wi,o,o + X] Vofl(^' 



5 '^0,0,-1 — ""-1,0 



"^0,0^ • 



i=l 



i=l 



into the initial symmetry. The expansion of this symmetry in a Taylor series gives an 
infinite number of symmetries. 

Application of the Gardner method for generation of conservation laws for the 
asymmetric equation 



a^io,oWi,o - /3(mo,oWo,i + wi,oMi, 



0. 



(41) 



is also considered. In this case we can again say that all necessary information for 
application of the method follows from multidimensional consistency. Equation fHTjl 
is consistent with the symmetric equation 



0. 



(42) 



We obtain a three dimensional consistent system by imposing fH2]) on 2 opposite sides 
of a cube and ( HTj) on the remaining 4 sides. From this consistency follows the Lax 
pair for ( 14T]) 



1 



M 



1 



/3mo,i 




Uo,o 
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Equation (HTl) satisfies all the criteria of integrability [22], namely there exist infinite 
number of conservation laws and symmetries, and a Lax pair. As is seen from this 
article all these properties are connected with multidimensional consistency. So, mul- 
tidimensional consistency is a universal criterion of integrability. It is not necessary 
for an equation to be consistent with itself, as we showed for ( HTl) it can be consistent 
with other equations. 

Here are some interesting topics for future research: 

• The conditions for the application of the Gardner method to the consistent equa- 
tion (equation and its BT). 

• The classification of quad-graph equations according to these conditions. 
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